arXiv:1503.07367vl [hep-ph] 25 Mar 2015 


Dilepton constraints in the Inert Doublet Model from Run 1 of the LHC 


Genevieve Belanger^, Beranger Dumont^’^, Andreas Goudelis^, Bjorn Herrmann^, Sabine Kraml^, Dipan Sengupta^ 
^LAPTh, Universite Savoie Mont Blanc, CNRS, B.P. 110, F-74941 Annecy-le-Vieux, France 
^Center for Theoretical Physics of the Universe, 

Institute for Basic Science (IBS), Daejeon 305-811, Republic of Korea 
^Laboratoire de Physique Subatomique et de Cosmologie, Universite Grenoble-Alpes, 

CNRS/IN2P3, 53 Avenue des Martyrs, F-38026 Grenoble, France 
'^Institute of High Energy Physics, Austrian Academy of Sciences, Nikolsdorfergasse 18, 1050 Vienna, Austria 

Searches in final states with two leptons plus missing transverse energy, targeting supersymmetric 
particles or invisible decays of the Higgs boson, were performed during Run 1 of the LHC. Recasting 
the results of these analyses in the context of the Inert Doublet Model (IDM) using Mad Analysis 5, 
we show that they provide constraints on inert scalars that significantly extend previous limits from 
LEP. Moreover, these LHC constraints allow to test the IDM in the limit of very small Higgs-inert 
scalar coupling, where the constraints from direct detection of dark matter and the invisible Higgs 
width vanish. 

PACS numbers: 12.60.-i, 12.60.Fr 


The Inert Doublet Model is one of the simplest ex¬ 
tensions of the Standard Model (SM). First introduced 
as a toy model for electroweak (EW) symmetry break¬ 
ing studies [1], it was later found to potentially improve 
naturalness [5], be compatible with Goleman-Weinberg 
EW symmetry breaking (EWSB) [3] as well as accom¬ 
modate dark matter (DM) through different potential 
mechanisms mm] and with interesting signatures 121- 
[U. It can moreover explain EW baryogenesis [iiiia, 
neutrino masses can be constrained from LEP mea¬ 
surements |15j and predicts a rich phenomenology at the 
LHG [ailMI]- 

The discovery of a Higgs-like particle with a mass of 
about 125 GeV at the LHG [5^ [33] has already made a 
tremendous impact on the phenomenology of the IDM, 
dramatically reducing its available parameter space [33] . 
Besides its interplay with Higgs phenomenology, the IDM 
can yield interesting signals at the LHG, involving in par¬ 
ticular two three [18] or multiple [19] leptons along 
with missing transverse energy, Such signals have 

so far been studied in the literature only as predictions 
and never as constraints from existing LHC results. In 
this work we show that existing data on -\- A™'®® 

searches from the 8 TeV LHC run — performed by the ex¬ 
perimental collaborations with supersymmetry (SUSY) 
or invisible Higgs decays in mind — begin to provide 
significant constraints on the IDM parameter space that 
are highly complementary to those obtained from DM 
observables. 

Let us begin by briefly presenting the IDM and set¬ 
ting some useful notations for the subsequent analysis. 
In the IDM, the SM is extended by the addition of a sec¬ 
ond scalar, $, transforming as a doublet under SU{2)l- 
This doublet $ is odd under a new discrete Z2 symmetry, 
whereas all other fields are even. In Feynman gauge we 
can write the two scalar doublets as 

+iG»)). *=(i(iA + a»)). 

( 1 ) 


where v = (0|iL|0) « 246 GeV denotes the vacuum 

expectation value of the neutral component of H. The 
h state corresponds to the physical SM-like Higgs boson, 
whereas and are the Goldstone bosons. The “in¬ 
ert” sector consists of a neutral CP-even scalar a 
pseudo-scalar A°, and a pair of charged scalars . The 
scalar potential of the model reads 


Uo = + A2|$|" 

+ + A4|i?^$P + y + h.c. 


( 2 ) 


The masses and interactions of the scalar sector are fixed 
by the scalar-potential parameters A1...5 and ^ 2 , which 
can be traded for the physically more intuitive set 


{mh, 

ruHO, 

rriAo, rnH±, Xl, 

A 2 }, (3) 

where the Higgs 

and inert scalar masses 

are given by 


ml 

— Ri -\- 3Aiu^, 

(4) 



= ^2 + Alu^, 

(5) 


m\o 

= rI + Asu^ 

(6) 



= l4 + 

(7) 

and the couplings Xl,s 

are defined as 



Al,s = 2 ("*^3 + A4 ± As). (8) 

The parameter is eliminated by using (after EWSB) 
the scalar potential minimization relation = 

2Aiu^. 

The Z2 symmetry forbids mixing among the compo¬ 
nents of H and $ and renders the lightest Z2-odd par¬ 
ticle stable which, if electrically neutral {i.e. or A°), 
can play the role of a DM candidate. Detailed accounts 
of the DM phenomenology of the IDM can be found, for 
example, in W- Sticking, for simplicity (but without 
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loss of generality), to the mass hierarchy < rriAo , 
the IDM can reproduce the observed [25] DM abundance 
according to the freeze-out mechanism in three regimes: 
the low-mass regime {rriffo < mw), where the relic den¬ 
sity is governed by the coupling Ai, of two H'^ particles 
to a Higgs boson and the distance of mno from mhl‘2- 
(the exact difference between and also plays 

a role when the former is larger than ~ 70 GeV), the 
intermediate-mass region (mw < mno ^ 115 GeV), 
where the relic density depends on m^o and Al, and 
the high-mass regime (which will not be of interest for 
this work) where all parameters of the scalar potential 
except A2 drastically affect the DM relic abundance. 

The IDM parameter space is subject to numerous con¬ 
straints that we shall impose throughout our analysis 
(for a detailed account see, e.g., [23]). First, there is a 
minimal set of theoretical requirements (stability of the 
EW vacuum, perturbativity of all couplings and pertur¬ 
bative unitarity of the scattering matrix) rendering the 
IDM a consistent and calculable quantum field theory 
(see also [26l[27]). We demand that these three require¬ 
ments be satisfied up to at least a scale of 10 TeV. We 
moreover have to ensure that there are no excessive con¬ 
tributions to the oblique parameters S', T and [/, for 
which we consider the 3a ranges from [28]. The analy¬ 
sis presented in m showed that (assuming for concrete¬ 
ness mfjo < m^a) the neutralino searches at LEP impose 
the bound ^ 100 GeV. A recast of the chargino 
searches at LEP [29] also amounts to a constraint on the 
charged scalar mass, namely mH± ^ mw- Finally, for 
m^o < mhl2, we have to make sure that the Higgs boson 
decay modes predicted by the IDM are in agreement to 
those observed at the LHC. For SM-like couplings to all 
SM particles, the invisible branching ratio of a Higgs into 
two particles is restricted to BR(/i —)■ inv.) < 0.12 at 
95% confidence level (CL) [30] (see also [3TIj34] L 

Additional constraints arise once the IDM is consid¬ 
ered as a DM model. Essentially the entire region where 
m^fo < mh/2 is ruled out due to the invisible Higgs 
branching ratio constraint, which imposes Al 6 x 10“^, 
while XENONIOO [3S] already eliminated the entire low- 
to-intermediate mass regime where the IDM can explain 
DM according to the thermal freeze-out mechanism. The 
only exception is a narrow region around tojjo ~ mu/2, 
which is characterised by A^ ~ 0, implying a feeble DM- 
Higgs (and, hence, DM-nucleon) coupling Pl]. The more 
recent LUX results [35] only render these constraints 
quantitatively stronger, forcing one to move even deeper 
into the “Higgs funnel” in order to find parameter space 
points compatible with the relic density constraints from 
Planck. 

All constraints coming either from invisible Higgs de¬ 
cays or from direct detection experiments however vanish 
in the limit Al —0. Therefore this regime of feeble Xl 
coupling appears to be extremely challenging experimen¬ 
tally. While it might be argued that, with the exception 
of the finely-tuned Higgs funnel region, this regime leads 
to a DM overabundance, it is definitely relevant to seek 


for constraints on DM models that do not depend in any 
way upon astrophysical or cosmological assumptions. In 
particular, the “no-go” argument for DM overabundance 
holds only in the context of a “vanilla” picture for the 
thermal history of the universe — numerous situations 
can arise in embedings of the IDM which can lead to 
an eventual dilution of the DM density, along the lines 
discussed, e.g., in m- Besides, the IDM needs not be 
viewed as a DM model, in the sense that it provides an 
interesting phenomenology in itself. In this spirit, we will 
see that the dilepton searches at the LHG provide con¬ 
straints that are complementary to those obtained from 
other experimental searches. 

The dilepton -|- signature in the IDM was first 

studied in m- Four processes provide the main signal 
contributions to this final state, namely 


qq 




It- 
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(9) 
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Process 0 directly depends on the coupling Al while 
all others involve only gauge couplings and thus depend 
only on the masses of the inert doublet fields. 

Searches for two opposite-sign leptons plus have 

been performed by both the ATLAS and CMS collabora¬ 
tions at Run 1. While no interpretation in the IDM was 
given by the experiments, there are various Higgs and 
SUSY analyses that are potentially sensitive to the IDM 
signatures. In the following, we will discuss the results 
obtained by ATLAS. 

On the SUSY side, there is a search for neutrali- 
nos, charginos, and sleptons [35]. In this context, the 
£~^£~ + signature arises from chargino-pair produc¬ 
tion followed by —>• or £^vlv£^ de¬ 

cays, or slepton-pair production followed by £^ —)■ 
decays. All signal regions (SRs) based on purely lep- 
tonic final states require |m« — mz\ > 10 GeV, i.e. they 
veto leptons coming from Z decays. One of the simpli¬ 
fied models considered in this search, x~^x~ VF+(—>■ 

£~^n)xiW~ (—>■ £~i')xi, can be directly matched with pro¬ 
cess ( [To] ). Several SRs are designed to optimize sensi¬ 
tivity to different mass splittings mj^± — m^o. Inter¬ 
estingly, note that the SUSY equivalent of process 
X2X1 — > > £^£ )xiXi, has not been considered in 

any ATLAS or CMS search. 

On the Higgs side, there is the search for invisible de¬ 
cays of a Higgs boson produced in association with a Z 
boson [35]. It requires |m« — mz\ <15 GeV, and can 
be matched to processes (111 and (12); it may also have 
sensitivity to process if m^o — m^jo > mz. Finally, 
there is another search in the £~^£~ final state, fo¬ 

cused on DM [30] ■ However, it requires > 150 GeV, 
which completely removes all of our signal. 
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In order to work out the current LHC constraints on 
the IDM, we recast these two ATLAS analyses using 
the MadAnalysis 5 [HI |H] framework. The SUSY 
search [35], was already available in the Public Analy¬ 
sis Database [13] as the recast code [H]. The invisi¬ 
ble Higgs search [39] was implemented and validated for 
this Letter and is available at |3S]. The signal gener¬ 
ation is done with MadGraph 5 m EZ] with model 
files generated using the FeynRules IDM implementa¬ 
tion presented in [24]; the particle widths are calculated 
with CalcHEP EHl HH]- The parton-level events are 
passed through Pythia 6.4 [SD] for parton showering and 
hadronization before being processed with the ‘MAStune’ 
version of Delphes 3 m (see Section 2.2 of [43]) for the 
simulation of detector effects. The number of events af¬ 
ter cuts are then evaluated with the recast codes jHJ US] . 
For the statistical interpretation, we make use of the 
module exclusion_CLs .py [43j : given the number of sig¬ 
nal, observed and expected background events, together 
with the background uncertainty, exclusion_CLs.py de¬ 
termines the most sensitive SR, the exclusion confidence 
level using the CLg prescription, and the nominal cross 
section (T 95 that is excluded at 95% Cljj 

The IDM parameter space, see Eq. (IS^is sampled tak¬ 
ing into account the following considerations. First, in 
light of the constraints discussed above which require A l 
to be tiny, process ( 11 ) is essentially irrelevant for the 
entire analysis. We can therefore choose A^, = 0 without 
loss of generality. Besides, A 2 is irrelevant for all observ¬ 
ables at tree-level. The mass of the charged inert scalar 
is important mostly for process ( 10 ), which comes with 
the price of an additional EW coupling factor with re¬ 
spect to (|^ and turns out to be numerically insignificant 
unless mH± is very light We are thus left with m^o and 
niHO to scan over. For mH ±, we choose two representa¬ 
tive values: TnH± = 85 GeV, which is the lower allowed 
limit by LEP, and mfi± =150 GeV, which is significantly 
higher but still safely within the bounds imposed by the 
T parameter, which limits the mass splitting between the 
inert scalar states (see also the analysis in [24 ] ). 

The main results of our analysis are presented in Fig¬ 
ure where we show fj, = (795 /(Tidm in the form of tem¬ 
perature plots in the (m^o, mno) plane for the two chosen 
values of mH± ■ Here, itidm is the cross section predicted 
by the model while CT 95 is the cross section excluded at 
95% CL. With this definition, regions where /i < 1 are 
excluded at 95% CL. 

As can be seen, the Run 1 ATLAS dilepton searches 
exclude, at 95% CL, inert scalar masses up to about 
35 GeV for pseudoscalar masses around 100 GeV, with 


^ Note that we do not simulate the backgrounds but take the back- 
ground numbers and uncertainties directly from the experimental 

publications m and m- 

^ Process in is also subdominant because the ZZH^H^ coupling 
is quadratic in the weak coupling whereas the coupling ZA^H^ 
is only linear. 
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FIG. 1. The ratio /r = o-gs/ciDM in the (m^o,m£fo) plane 
for two representative values of the charged inert scalar mass, 
'mu± = 85 GeV (upper panel) and = 150 GeV (lower 

panel). The solid black lines are the 95% GL exclusion con¬ 
tours, fj, — 1. The dashed black lines are given for illustration 
and correspond to the fx = 0.5 and = 2 contours. The grey 
dashed lines indicate m^o — m^jo = mz- 


the limits becoming stronger for larger m^o, reaching 
« 45 (55) GeV for ~ 140 (145) GeV and mH± = 85 
(150) GeV. For massless A° masses up to about 135- 
140 GeV are excluded (note that uiho and ttiao are gen¬ 
erally interchangeable here). Several interesting features 
merit some discussion. 

First, we observe that the constraints are slightly 
stronger for heavier charged scalars. This is in part due to 
the small contribution from process ( 10 ) and from qq —>■ 
—?► AH^ —H where one of the leptons 
is missed: although the cross section is much larger for 
mH± = 85 GeV as compared to mH± = 150 GeV, the 
resulting leptons are much softer and almost never pass 
the signal requirements. A more significant difference be¬ 
tween the mH± = 85 GeV and 150 GeV cases arises from 
the fact that at large m^o the signal from process ([^ 
is suppressed by the decay A° —>• followed by 

—>■ W^H^, that competes with A° —>■ While 

the former decay mode also leads to dileptons, these lep¬ 
tons are, as above, much softer and almost never pass 
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m^± = 150 GeV 



m^o (GeV) 


FIG. 2. Most sensitive SR in the {m^o,m,uo) plane. Also 
shown as solid lines are the individual 95% CL exclusion 
curves obtained from the SUSY (in red) and Higgs (in or¬ 
ange) analyses considered in this study. 

the signal requirements. For these reasons, the rnjj± de¬ 
pendence is well captured by Figure increasing mff± 
further to, e.g., 300 GeV has very little effect]^ 

Second, and more importantly, the limits on rriHO be¬ 
come stronger for larger masses (above the diago¬ 
nal line of rrij^o = -|- mz)- This is also a kine- 

matical effect, because the leptons originating from the 
—)■ decay become harder when there is more 

phase space available and thus pass the signal selection 
cuts more easily. Similarly, for small mass splittings be¬ 
tween iL° and A°, the produced dileptons are softer for 
smaller m^o. 

Third, once the mass splitting (m^o — rriHo) is large 
enough for the Z to be on-shell, the Z veto in the SUSY 
analysis eliminates most of the signal. In this region it 
is the Zh —>■ -\- if™*®® search that gives the stronger 

limit. This is illustrated in Figure which shows the 
most sensitive SRs in the (m^o, rrij^o ) plane. The con¬ 
crete example used is for = 150 GeV but the results 
basically do not depend on the mass of the charged scalar. 
As expected, for — rriffo < mz the SUSY dilepton 
search is more sensitive, while for m^o — m^o > mz ex¬ 
clusion is driven by the Zh, h —>■ inv. search. Also shown 
are the individual 95% CL exclusion curves obtained from 
the SUSY and Higgs analyses considered in this study. It 
is worth noting that there is a small overlap region at low 
mass, m^o < 25 GeV, where the SUSY search also 
gives a limit despite the Z boson from the decay being 
on-shell. This is due to a sizeable tail of the mu distribu¬ 
tion extending below mu = mz — 10 GeV that allows to 
pick up some signal events in the SUSY SR WWb (this 


® In the strict IDM framework, mjj± = 300 GeV is excluded by 
EW precision tests for the {m^o region depicted in Fig- 


tail was already noticed in [H])- 

The evolution of these constraints in the next run of the 
LHC is a non-trivial question. A naive rescaling of signal 
and background numbers (see, e.g., [53]), assuming that 
the acceptanceX efficiency values remain the same, gives 
that at 13 TeV and an integrated luminosity of 100 fb“^, 
the 95% CL reach should extend up to ^ « 1.2 (1.6) 
above (below) the dashed grey line in Figure This im¬ 
proves to /i « 2.1 (2.7) with 300 fb“^, thus covering a 
major part of the (TO^o,m//o) plane shown in Figure]^ 
The exact reach will depend on many unknown factors, 
including the performance of the LHC accelerator and of 
the detectors at Run 2, improvements in the identifica¬ 
tion and reconstruction algorithms, a better estimation 
of the SM background, and modified selection cuts. The 
important point we want to make is that even with anal¬ 
yses that are not directly targeted at the IDM, the cov¬ 
erage of the parameter space is quite good. On the other 
hand, the LHC sensitivity could certainly be improved 
with an analysis optimized for pp —>• A°H° —^ Z^*'>H^H^ 
(or alternatively pp —>■ X 2 X 1 In particular, 

angular separation variables as discussed in m should 
prove useful to enhance the sensitivity. We therefore 
encourage the experimental collaborations to perform a 
dedicated search for inert scalars at 13 TeV. 

In conclusion, the results from Run 1 of the LHC pro¬ 
vide relevant constraints on the IDM model that signif¬ 
icantly extend previous limits from LEP. They are com¬ 
plementary to the cosmological and astrophysical con¬ 
straints, as they allow to test the model in the limit 
Al —)■ 0 , where the constraints from the direct detec¬ 
tion of DM and the invisible Higgs width vanish. It is 
also worth noting that these limits are independent of 
DM considerations, that may be different for extensions 
of the model like in l5H - f57] . The 95% CL limits we 
derived from the dilepton -|- U™'®® Higgs and SUSY anal¬ 
yses exclude inert scalar masses of up to about 55 GeV in 
the best cases. It should be possible to push these con¬ 
straints beyond mt/^ « 62.5 GeV at Run 2 of the LHC, 
thus starting to probe the Higgs funnel region of the IDM, 
a regime that is very hard to test through direct DM 
searches. We stress that these conclusions were reached 
based on existing analyses that were not optimised at all 
for the IDM signal. A dedicated experimental analysis 
for the IDM at 13 TeV would be highly desirable. 

Acknowledgements: BD and SK thank Maria 
Krawczyk and Dorota Sokolowska for useful discussions. 
This work was supported in part by the ANR project 
DMAstroLHC, ANR-12-BS05-0006, the “Investisse- 
ments d’avenir, Labex ENIGMASS”, and by the IBS un¬ 
der Project Code IBS-R0I8-DI. AG is supported by the 
New Frontiers program of the Austrian Academy of Sci¬ 
ences. 


urefTl but can be allowed in extensions of the model like in m 













5 


[1] N. G. Deshpande and E. Ma, Phys. Rev. D 18, 2574 
(1978) 

[2] R. Barbieri, L. J. Hall, and V. S. Rychkov, Phys. Rev. D 
74, 015007 (2006), arXiv:hep-ph/0603188 [hep-ph] 

[3] T. Hambye and M. H. Tytgat, Phys. Lett. B 659, 651 
(2008), arXiv:0707.0633 [hep-ph] 

[4] L. Lopez Honorez, E. Nezri, J. E. Oliver, and M. H. 
Tytgat, JCAP 0702, 028 (2007), arXiv:hep-ph/0612275 
[hep-ph] 

[5] T. Hambye, F.-S. Ling, L. Lopez Honorez, and J. Rocher, 
JHEP 0907, 090 (2009), arXiv:0903.4010 [hep-ph] 

[6] L. Lopez Honorez and C. E. Yaguna, JCAP 1101, 002 
(2011), arXiv:1011.1411 [hep-ph] 

[7] M. Gustafsson, E. Lundstrom, L. Bergstrom, and 
J. Edsjo, Phys. Rev. Lett. 99, 041301 (2007), arXiv:astro- 
ph/0703512 [astro-ph] 

[8] P. Agrawal, E. M. Dolle, and C. A. Krenke, Phys. Rev. D 
79, 015015 (2009), arXiv:0811.1798 [hep-ph] 

[9] S. Andreas, M. H. Tytgat, and Q. Swillens, JCAP 0904, 
004 (2009), arXiv:0901.1750 [hep-ph] 

[10] E. Nezri, M. H. Tytgat, and G. Vertongen, JCAP 0904, 
014 (2009), arXiv:0901.2556 [hep-ph] 

[11] C. Garcia-Cely and A. Ibarra, JCAP 1309, 025 (2013) 
arXiv:1306.4681 [hep-ph] 

[12] T. A. Chowdhury, M. Nemevsek, G. Senjanovic, and 
Y. Zhang, JCAP 1202, 029 (2012), arXiv: 1110.5334 [hep- 
ph] 

[13] D. Borah and J. M. Cline, Phys.Rev. D86, 055001 
(2012), arXiv: 1204.4722 [hep-ph] 

[14] E. Ma, Phys. Rev. D 73, 077301 (2006), arXiv:hep- 
ph/0601225 [hep-ph] 

[15] E. Lundstrom, M. Gustafsson, and J. Edsjo, 
Phys. Rev. D 79, 035013 (2009), arXiv:0810.3924 [hep- 
ph] 

[16] Q.-H. Cao, E. Ma, and G. Rajasekaran, Phys. Rev. D 
76, 095011 (2007), arXiv:0708.2939 [hep-ph] 

[17] E. Dolle, X. Miao, S. Su, and B. Thomas, Phys. Rev. D 
81, 035003 (2010), arXiv:0909.3094 [hep-ph] 

[18] X. Miao, S. Su, and B. Thomas, Phys. Rev. D 82, 035009 
(2010), arXiv: 1005.0090 [hep-ph] 

[19] M. Gustafsson, S. Rydbeck, L. Lopez-Honorez, 
and E. Lundstrom, Phys.Rev. D86, 075019 (2012) 
arXiv:1206.6316 [hep-ph] 

[20] A. Arhrib, R. Benbrik, and N. Gaur, Phys. Rev. D 85, 


095021 (2012), 

arXiv:1201.2644 [hep-ph] 

[21] B. Swiezewska and M. Krawczyk, Phys.Rev. D88, 035019 

(2013), 

arXiv:1212.4100 [hep-ph] 


[22] G. Aad et al. (ATLAS Collaboration), Phys. Lett. B 716, 
1 (2012), arXiv:1207.7214 [hep-ex] 


[23] S. Chatrchyan et al. (CMS Collaboration), Phys. Lett. B 
716, 30 (2012), arXiv:1207.7235 [hep-ex] 

[24] A. Goudelis, B. Herrmann, and O. Stal, JHEP 1309, 
106 (2013), arXiv: 1303.3010 [hep-ph] 

[25] P. Ade et al. (Planck), Astron.Astrophys. 571, A16 
(2014), arXiv: 1303.5076 [astro-ph.CO] 

[26] N. Khan and S. Rakshit, (2015), arXiv:1503.03085 [hep- 
ph] 

[27] B. Swiezewska, (2015), arXiv: 1503.07078 [hep-ph] 

[28] M. Baak, M. Goebel, J. Haller, A. Hoecker, D. Ludwig, 
et al, Eur.Phys.J. C72, 2003 (2012), arXiv: 1107.0975 
[hep-ph] 


[29] A. Pierce and J. Thaler, JHEP 0708, 026 (2007) 
arXiv:hep-ph/0703056 [HEP-PH] 

[30] J. Bernon, B. Dumont, and S. Kraml, Phys.Rev. D90, 
071301 (2014), arXiv:1409.1588 [hep-ph] 

]31] G. Belanger, B. Dumont, U. Ellwanger, J. Gunion, and 

S. Kraml, Phys.Lett. B723, 340 (2013), arXiv: 1302.5694 
[hep-ph] 

]32] A. Falkowski, F. Riva, and A. Urbano, JHEP 1311, 111 
(2013), arXiv:1303.1812 [hep-ph] 

]33] G. Belanger, B. Dumont, U. Ellwanger, J. Gu¬ 
nion, and S. Kraml, Phys.Rev. D88, 075008 (2013) 
arXiv: 1306.2941 [hep-ph] 

[34] P. Bechtle, S. Heinemeyer, O. Stal, T. Stefaniak, and 
G. Weiglein, JHEP 1411, 039 (2014), arXiv: 1403.1582 
[hep-ph] 

[35] E. Aprile et al. (XENONIOO Collaboration), 
Phys.Rev.Lett. 109, 181301 (2012), arXiv: 1207.5988 
[astro-ph. CO] 

[36] D. Akerib et al. (LUX Collaboration), Phys.Rev.Lett. 
112, 091303 (2014), arXiv:1310.8214 [astro-ph.CO] 

[37] G. Gelmini, P. Gondolo, A. Soldatenko, and C. E. 
Yaguna, Phys.Rev. D74, 083514 (2006), arXiv:hep- 
ph/0605016 [hep-ph] 

[38] G. Aad et al. (ATLAS Collaboration), JHEP 1405, 071 
(2014), arXiv:1403.5294 [hep-ex] 

[39] G. Aad et al. (ATLAS Collaboration), Phys.Rev.Lett. 
112, 201802 (2014), arXiv: 1402.3244 [hep-ex] 

[40] G. Aad et al. (ATLAS Collaboration), Phys.Rev. D90, 
012004 (2014), arXiv:1404.0051 [hep-ex] 

]41] E. Conte, B. Fuks, and G. Serret, Corn- 
put.Phys.Commun. 184, 222 (2013), arXiv: 1206.1599 
[hep-ph] 

[42] E. Conte, B. Dumont, B. Fuks, and C. Wymant, 
Eur.Phys.J. C74, 3103 (2014), arXiv: 1405.3982 [hep-ph] 

[43] B. Dumont, B. Fuks, S. Kraml, et al., Eur.Phys.J. C75, 
56 (2015), arXiv: 1407.3278 [hep-ph] 

]44] B. Dumont, (2014), 10.7484/INSPIRE- 

HEP. DATA. HLMR.T56W. 2 

[45] B. Dumont, (2015), 10.7484/INSPIRE- 

HEP. DATA. RT3V.9PJK 

[46] J. Alwall, M. Herquet, F. Maltoni, O. Mattelaer, and 

T. Stelzer, JHEP 1106, 128 (2011), arXiv: 1106.0522 
[hep-ph] 

]47] J. Alwall, R. Frederix, S. Frixione, V. Hirschi, F. Maltoni, 
et al., (2014), arXiv:1405.0301 ]hep-ph] 


[48] A. Pukhov, (2004), arXiv:hep-ph/0412191 ]hep-ph] 


[49] A. Belyaev, N. D. Christensen, and A 

. Pukhov, Com- 

put.Phys.Commun. 184, 1729 (2013) 

arXiv: 1207.6082 


[hep-ph] 

[50] T. Sjostrand, S. Mrenna, and P. Z. Skands, JHEP 0605, 
026 (2006), arXiv:hep-ph/0603175 [hep-ph] 

[51] J. de Favereau et al. (DELPHES 3), JHEP 1402, 057 
(2014), arXiv:1307.6346 [hep-ex] 

]52] G. Belanger, K. Kannike, A. Pukhov, and M. Raidal, 
JCAP 1406, 021 (2014), arXiv: 1403.4960 [hep-ph] 

[53] A. Bharucha, S. Heinemeyer, and F. von der Pahlen, 
Eur.Phys.J. C73, 2629 (2013), arXiv:1307.4237, 

[54] A. D. Banik and D. Majumdar, Eur.Phys.J. C74, 3142 
(2014), arXiv:1404.5840 [hep-ph] 

[55] P. Ko, Y. Omura, and C. Yu, JHEP 1411, 054 (2014) 
arXiv:1405.2138 [hep-ph] 




6 


[56] M. Hindmarsh, R. Kirk, J. M. No, and S. M. West, [57] C. Bonilla, D. Sokolowska, J. L. Diaz-Cruz, M. Krawczyk, 
(2014), arXiv:1412.4821 [hep-ph] and N. Darvishi, (2014), arXiv: 1412.8730 [hep-ph] 



